version 2 



Dynamics of Air-Fluidized Granular System Measured by the 

Modulated Gradient Spin-echo 



Janez Stepisnik/'O Samo Lasic, 1 Igor Sersa, 2 Ales Mohoric, 1 and Gorazd Planinsic 1 

1 University of Ljubljana, Faculty of Mathematics and Physics, 
Physics Department, Jadranska 19, 1000 Ljubljana, Slovenia^ 
2 Josef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia 
(Dated: February 2, 2008) 

Abstract 

The power spectrum of displacement fluctuation of beads in the air-fluidized granular system is 
measured by a novel NMR technique of modulated gradient spin-echo. The results of measurement 
together with the related spectrum of the velocity fluctuation autocorrelation function fit well to 
an empiric formula based on to the model of bead caging between nearest neighbours; the cage 
breaks up after a few collisions The fit yields the characteristic collision time, the size of bead 
caging and the diffusion-like constant for different degrees of system fluidization. The resulting 
mean squared displacement increases proportionally to the second power of time in the short- 
time ballistic regime and increases linearly with time in the long-time diffusion regime as already 
confirmed by other experiments and simulations. 

PACS numbers: 45.70.Mg, 76.60.Lz 
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INTRODUCTION 



Sand dunes, grain silos, building materials, catalytic beds, filtration towers, riverbeds, 
snowfields, and many foods are granular systems consisting of large number of randomly 
arranged macroscopic grains. Despite their apparent simplicity granular materials exhibit 
a host of unusual behaviours, whose unravelling more often than not appears to challenge 
existing wisdom of science 

Fluidized granular bed is a system of randomly arranged, macroscopic grains in which 
the driving force of motion is container shaking or gas flow through the granular system. Al- 
though, these systems are of tremendous technological importance in catalysis of gas-phase 
reactions, transport of powders, combustion of ores, and several other industrial processes, 
we do not have sufficient understanding of the fluid state of granular medium that is analo- 
gous to macroscopic properties of liquids. Two particularly important aspects contribute to 
the unique properties of granular materials: thermodynamics plays no role, and interactions 
between the grains are dissipative, because of static friction and inelasticity of collisions. 
Several theoretical efforts start towards building granular fluid mechanics by considering 
the medium as a dense, inelastic gas with the temperature defined by induced local velocity 
fluctuations 0, 0] . The autocorrelation function of velocity fluctuation is the basis of many 
thermodynamic and hydrodynamic models which aim to provide a statistical description of 
a granular system in terms of a single particle dynamic. In the simulation of hard-sphere 
fluid, Alder and Wainwright f| found a strong dependence of the velocity autocorrelation 
function (VAF) on the system density. Only for very low particle densities, it decays expo- 
nentially with the Enskog correlation time jj], while at higher densities, a negative long-time 
tail appears as the result of caging by adjacent spheres. 

Although, the experimental techniques used to study the motion of granular systems span 
a wide range of approaches and sophistication, very few attempts were made to finger into 
the details of grain motion. To the best of our knowledge, only the tracking of bead motion 
by the positron emission [q] and by the CCD camera iOj were able to glimpse to the VAF of 
a fluidized granular bed. 

The NMR gradient spin-echo is a tool that yields not only macroscopic but also mi- 



croscopic dynamic variab 
of spin bearing particles 



es due to the relation between spin-echo attenuation and VAF 
lfll fllL \v\ . However, this potential of spin-echo has been only 
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particularly, when used for the study of granular 



partially exploited 
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In this letter, we report on the first application of modulated gradient spin-echo (MGSE) 
to measure autocorrelation spectra of motion in a fluidized granular bed. MGSE is a method, 
in which a repetitive train of radiofrequency (RF) pulses with interspersed magnetic field 
gradient pulses or a gradient waveform periodically modulates the spin phase, so that the 
spin-echo attenuation is proportional to the power spectrum of the spin displacement fluc- 
tuation (PSDF) or the spectrum of VAF. The frequency range of the measurement is de- 
termined by the rate of spin phase modulation, which can be between a few Hz to about 
a few 10 kHz at the present state of NMR hardware. This covers the expected range of 
fluidized bead motion Q|. Generally, the measurement of the fluctuation autocorrelation 
by the MGSE method can be considered as a low-frequency complement to the non-elastic 
neutron scattering method that covers the range above GHz frequency. 




FIG. 1: The high-speed camera snap-shot of grains moving in the air fluidized system. The curves 
show the paths of two representative grains during 20 ms long MGSE sequence. The arrows at the 
bottom indicate the air flow inlet through the perforated hose. 



When using NMR spin-echo to detect the translational motion of a spin-bearing particle, 
we observe the cumulative effect of small perturbations of spin precession frequency due to 
the displacements of the spin-bearing particle in the magnetic field gradient. The random 
collisions between grains in the fluidized granular system shift a bead from the path of mean 
motion driven, for example, by the flow of air. As the phase of the spin-echo signal depends 
on the displacements along the mean path of motion, the signal attenuation depends on the 
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PSDF, I z (oo) , as [11 

(3{t) = 1t h^)\G{u,r)\ 2 d^ (1) 

whenever the displacement fluctuations from the mean path are short compared to the phase 
shift grating created by the applied gradient [2JJ . The spectrum of the effective gradient 
G e ff(t) is defined as ll| 

G(co,r)= f G eff (t)e-^dt, (2) 

<j 

where r is the time of spin-echo measured from the first excitation pulse. The effective 
gradient changes the sign upon every application of the 7T-RF pulse. According to the 
Wiener-Khintchine, the PSDF is related to the displacement autocorrelation function j^] 

I z ( u ) = - r (Az(t) Az(0)}e- lut dt, (3) 

7T J-oo 

where Az(t) = z{t) — (z(t)) denotes the displacement fluctuation along the direction of 
applied gradient with {z(t)) being the mean spin location at time t. The PSDF is related 
to the power spectrum of the velocity fluctuation or the spectrum of the VAF 

D(u) = I z (oo)oo 2 , (4) 

and to the mean squared displacement of grain fluctuation 

{[Az(t) - Az(0)] 2 } = - / I z (uj)(l-cos(ujt))duj. (5) 

n Jo 



The fluctuating part of displacement cannot be distinguished from the mean displace 
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ment with the use of the conventional spin-echo with two gradient pulses (PGSE) 
where the resulting propagator includes the average of a grain flow and grain displacement 
fluctuations. The flow compensating four gradient pulse variant of PGSE j3| is unsuitable 



for the PSDF sampling, because its gradient spectrum covers a broad-frequency range 24]. 
However, a RF-gradient sequence, that causes a cyclic spin phase modulation, yields a gra- 
dient spectrum with sharp peaks. After iV modulation cycles of the period T, the resulting 
spin-echo attenuation is 

(3{NT) = 2 1 2 NTY,L(nuj m )\g(nuj m )\ 2 . (6) 
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FIG. 2: RF- gradient sequence that yields the spectrum of effective gradient with the dominant 
peak at the frequency uj = 2n/T. 



Here I z (nuj m ) are values of the PSDF at the peaks of the gradient spectrum at ui = with 



i/ T G eff (t)e- iujt dt, while the peak 



n = 0, ±1, ±2, ±3.... The peak amplitudes are g(co 

widths are 2n/NT. A sequence with a single dominant peak in the gradient spectrum 
can be designed with a proper phase cycling. The CPMG train of it RF pulses |26;] applied 
to spins in a constant magnetic field gradient with magnitude G forms dominant peaks at 
uj m = ±2tt/T, whenever the train of 7T-RF pulses, which are repeated at the T/2 intervals, 
is applied a quarter of the period after the excitation 7r/2-RF pulse, as shown in Fig|21 The 
resulting spin-echo attenuation at r = NT 

87 2 G 2 NT 



P{NT,u r , 



7T^ 



(7) 



can be used to sample the PSDF by varying the modulation period T . 



EXPERIMENTAL PROCEDURE 



The experiment was carried out on a TecMag NMR spectrometer with a 2.35 T hori- 
zontal bore superconductive magnet. The spectrometer was equipped with micro-imaging 
accessories and with reversed Helmholtz gradient coils with 0.25 T/m peak magnetic field 
gradient. 

The described MGSE method was applied to a granular system made up of 100 pharma- 
ceutical 3-mm diameter, oil-filled, hard plastic spherical beads with the restitution coefficient 
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0.85 placed in a cylindrical chamber. The container, build from a piece of a plastic syringe 
tube of 26 mm length and 23.5 mm diameter, was placed inside the RF coil with the cylinder 
axis being horizontal. A perforated hose of 5 mm outer diameter, attached to the bottom 
inner wall of the container, is used as a diffuser of air up-flow as shown in FigO Several 
holes of 0.5 mm diameter are drilled uniformly across the surface of the container to serve as 
the air outlets. The degree of ffuidization was regulated by the air pressure in the hose and 
monitored by a high speed camera (600 fps) through the transparent walls of the cylindrical 
container. At applied air pressures, we observed that initial tendency of uplifted beads, to 
circulate around the container walls, breaks down by collisions and dispersed air flow lead- 
ing to a random motion of beads. The system was considered as a fully fluidized because 
a distinct bead clustering or bubbling was not observed (^tJ . FigQ shows two paths of 
representative single grains at the air pressure of 0.5 bar during 10 steps that correspond to 
20 ms. The paths of beads in the bulk are short enough that the collisions with the container 
walls are seldom during the measuring interval. 

The PSDF of fluidized system was measured in the static gradient field of 0.0625 T/m 
in the vertical direction, perpendicular to the axis of cylindrical container. The axis was 
parallel to the static magnetic field. The amplitudes of spin-echo were recorded at different 
modulation periods, T, while keeping the acquisition time constant, t = NT = 20 ms. The 
co-addition of echoes, taken at different steps of RF phase cycle of CPMG train, is used to 
remove the echoes of spurious coherence path-ways as well as a residual phase shift of the 
grain flow. Thus, the applied method was able to detect only the displacement fluctuation 
caused by the collisions and the air flow. 

RESULTS AND DISCUSSION 

The MGSE measurement gives I z (uj) of beads in the air-fluidized system at different gas 
pressures, as shown in FigJSJ The spectra exhibit a clear exponential decay at frequencies 
above 400 Hz, which is more distinctive for stronger air- flow at 0.5 bar. Evidently, faster 
grain motion increases the spin-echo attenuation at higher modulation frequencies magnify- 
ing the high-frequency range of measurement, but hindering the examination below 300 Hz, 
because of a weak signal. At the reduced air-pressure of 0.25 bar, the MGSE measurement 
can trace I z (u) below 400 Hz, yielding the frequency dependence that passes from an ex- 
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FIG. 3: The displacement fluctuation spectra of oil-filled beads fluidized by the air-flow at pressures 
of 0.25 bar (empty circles) and 0.5 bar(full circles) . At higher frequencies, the experimental points 
indicate clear exponential decay. At the low pressure, there is a very distinctive ^-dependence at 
low modulation frequencies. Curves display the best fit by the empiric formula (EqJHJ). 

ponential at high frequencies into a 1/u 2 low-frequency dependence. The power spectrum 
of VAF, D(uj) in Fig0J\, calculated accor ding to Eq0J is similar to the power spectrum of 
over-damped thermal harmonic oscillators [28| . The form of the spectrum is not Lorentzian 
as one would expect by assuming the Enskog exponential decay of the velocity autocorre- 
lation function. Interestingly, the positron emission measurements give a similar spectral 
lobe for the vibrofluidized granular bed |8|], which is explained as a bead caging within the 
experimental cell. Better signal to noise ratio of MGSE measurements enables a detailed 
analysis of I z {uj) and related D{uj). In Fig0]\l, D(u>) shows a clear u; 2 -dependence at low 
frequencies, which is typical for restricted diffusion. In restricted diffusion case, the slope 
of D{uj) vs. uj 2 plot provides information about the size of spin confinement, while the 
intersection with the ordinate, -D(O), gives the inter-pore diffusion rate in a system with 
inter-pore channels or permeable walls. Thus, a small value of -D(O) in Fig^A.1 could be 
interpreted by the model of bead ballistic motion between successive collisions, where many 
collisions are required to break out the caging formed by adjacent beads. The exponential 
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frequency dependence of I z (uj) in the high-frequency range, as observed in our experiments, 



29|. However, 



is different from the Lorentzian spectrum for diffusion in porous media 
the length of ballistic grain motion depends on the distribution of adjacent beads, which 
might be responsible for the observed form of I z (uj) at high frequencies. By taking into 
account the experimentally observed character of the spectrum, an empiric formula 

LO 

is introduced, which gives the low-frequency dependence similar to the restricted diffusion 
and the exponential decay at high frequencies. Here J (£ 2 ) is the size of the bead 



caging, t c is the characteristic time of ballistic motion and D is the diffusion-like constant 
describing hopping between cages. The mean squared displacement, calculated from this 
formula according to Eq|5J is (Az 2 (t)) « 4(£ 2 )t 2 /7rr 2 in the short-time ballistic reg ime and 
is (Az 2 (t)) ~ 2Dt in the long-time diffusion regime. The optical measurements 30] and the 
simulations [31] give a similar time- dependence. 

The curves in Fig 03 correspond to the best fits of the empiric formula to the experimental 
data. For the air-flow at 0.25 bar, the fitting parameters are: (£ 2 ) = 1.2 x 10 -9 m 2 , r c = 0.22 
ms and D = 0.47 x 10~ 6 m 2 /s, with the relative error of 5%. The air-flow at 0.5 bar gives 
a very clear exponential dependence with the fitting parameters: (£ 2 ) = 4.4 x 10~ 9 m 2 , and 
t c = 0.36 ms, but less exact D = 0.6 x 10 -6 m 2 /s, because of the attenuation cut-off in 
the low-frequencies range. The characteristic time of ballistic motion related to the bead 
collision frequency can be obtained from the position of D (u) -maximum at uo = 2/r c . 

The Fourier transform of D(u>), modelled with the above fitting parameters, gives the 
VAF of fluctuation, (Av z (t) Av z (0)), as shown in Fig0j3. Its intersection with the ordinate 
provides the mean squared of the velocity fluctuation, which is by definition proportional 
to the temperature of the fluidized system: T ~ (v 2 (0)). It demonstrates an increase of the 
system temperature by about 45% as the pressure, driving the air-flow, increases from 0.25 
bar to 0.5 bar. 



CONCLUSION 

The new MGSE method provides the PSDF and the spectrum of VAF of beads in the air 
fluidized granular system. The data fit well to the empiric formula that corresponds to the 
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FIG. 4: A.) The power spectrum of grain velocity fluctuation in the system of oil-filled beads 
fluidized by the air-flow at 0.25 bar (squares) and 0.5 bar (crosses). Al.) The intersection of line 
with the ordinate in D(uj) vs. w 2 -plot gives the diffusion-like constant. B.) The Fourier transform 
of D{uo) gives the VAF of beads for air flow at the pressure 0.25 bar (dashed) and 0.5 bar. 

model of bead caging within the space of adjacent colliding beads, which breaks up after 
many collisions. In references j^, , the positron emission measurement result in a similar 
spectra shape of VAF for the grain motion in a vibrofluidized granular system. Similar 
VAF, as shown in Fig0j3, was measured in the shear flow of a granular system by using 
a CCD camera e 9], where a break down of the exponential decay at high particle densities 
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was reported. In the statistical sense, the observed exponential form of the PSDF could 
be related to the observed exponential distribution of forces in a sheared granular bed as 
reported in Ref. Q|. 
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